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ABSTRACT 

The Hubble Ultra Deep Field Survey has not only provided the deepest optical 
and near infrared views of universe, but has enabled a search for the most distant 
supernovae to z ~ 2.2. We have found four supernovae by searching spans of 
integrations of the Ultra Deep Field and the Ultra Deep Field Parallels taken 
with the Hubble Space Telescope paired with the Advanced Gamera for Surveys 
and the Near Infrared Multi Object Spectrometer. Interestingly, none of these 
supernovae were at z > 1.4, despite the substantially increased sensitivity per 
unit area to such objects over the Great Observatories Origins Deep Survey. We 
present the optical photometric data for the four supernovae. We also show that 
the low frequency of Type la supernovae observed at z > 1.4 is statistically 
consistent with current estimates of the global star formation history combined 
with the non-trivial assembly time of SN la progenitors. 

Subject headings: Surveys — supernovae: general — supernovae:individual (K0302- 
001, SN 20031t, SN 20031u, SN 2004R 


1. Introduction 

The Hubble Ultra Deep Field (UDF) survey (GO 9978, S. Beckwith, PI) is the deepest 
view of the universe in the optical and near infrared obtained to date. Imaging was acquired 
with the Advanced Gamera for Surveys (AGS) at a position located within the Great Ob¬ 
servatories Origins Deep Survey (GOODS, Giavalisco et al. 2004a) South Field, centered at 
R. A. (J2000)=03h32m39’00 and Deck (J2000)=-27°4r29.1". A total of 400 RFTorbits were 


^ Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555 
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accumulated to achieve a lOa detection threshold of approximately 29 mag (Vega) in the 
F435hh, F6061V, and F775hh-bands, and approximately 28 mag in the F850LF-band. To 
further enhance the vast array of deep multi-wavelength imaging for this region, deep near 
infrared images were also obtained in the UDF target held (IRUDF; GO 9803, R. Thomp¬ 
son, PI) by acquiring individual pointings with the Near Infrared Camera and Multi Object 
Spectrometer (NICMOS) camera 3, each 8 orbits in depth in both the FllOhF and F160IF- 
bands, to cover a 3 x 3 mosaic near the center of the UDF. This observation was divided into 
two epochs separated by 90 days to enable searches for SNe. Another extremely deep optical 
held, an ACS “parallel held”, observed during the imaging of the IRUDF, overlapped with 
GOODS South observations. The UDF, the UDF ACS Parallel (UDFP), and the IRUDF 
surveys are unique in that they are the hrst surveys sensitive to Type la supernovae (SNe la) 
to z < 2.2. 

With an elasped time of less than 4 Gyr between the hrst generation of stars {z ~ 100) 
and SNe la at z > 1.4, the observed number of these events at 1.4 < 2 ; < 2.2 could provide an 
important probe of the assembly time required by SNe la progenitors (Strolger et ah 2004, 
hereafter S04). The rate in which SN la events occur is governed by the rate in which their 
progenitor stars form, and the time required for the SN la progenitor to develop into a SN 
la event. In the framework of cosmic time, these components are the star formation rate 
history, and the population assembly time, or the delay time function (Dahlen & Fransson 
1999, S04). Together, the star formation rate history and the delay time function describe 
a model for the SN la rate history [JZia{z)] , which can be compared to observations of the 
SN la rate in different redshift regimes. 

In principle, the observed discovery rate of SNe at z > lA could help to distinguish 
between viable models of TZja{z). In S04, it was shown that observations of the SN la 
rate over a large redshift range (encompassing the most of the last ~ 10 billion years of 
the universe) are consistent with the combination of the star formation rate history from 
rest-frame U-band galaxy studies [hereafter SFR(/( 2 ;), see Giavalisco et al. (2004b)], and a 
Gaussian distribution of delay times, with a mean ~ 4 Gyr and a dispersion of ~ 1 Gyr. 

However, recently it has been shown that the rate SNe la per galaxy at < 0.1 does 
appear to increase towards later galaxy types, e. g. the SN la rate in irregular (dwarf) galaxies 
is approximately 10 times larger than in elliptical galaxies of the same total mass (Mannucci 
et al. 2004). This would seem to imply that the SN la rate more closely traces the star 
formation rate history rather than being largely delayed from it. This is difficult to resolve 
in light of the GOODS supernova data. However, one could hypothesize that the rest-frame 
U-band galaxy observations (including dust corrections), in reality, provide a rather poor 
tracer of actual star formation history, or only a lower limit (Ghary & Elbaz 2001). Or, 
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perhaps the completeness thresholds for the GOODS SN survey, the sensitivity as a function 
of change in flux, were signihcantly overestimated. 

A simple test of the S04 model would be to push the limits of survey sensitivity well 
beyond those achieved in the GOODS survey, and continue to look for high redshift SNe la 
at z > 1.4. One can then compare observed yields to those expected in this redshift regime 
from the S04 and other viable models of the SN la rate history. If 7lia{z) is indeed delayed 
by 4 Gyr from the SFR(/( 0 ), then it is expected that the rate at z > 1.4 will be very low, 
and thus essentially no (<^ 1) SNe la in this range could be found in this deep survey. Even 
if the Tiia{z) is closely tied to the SFRt/( 2 :), i. e. with mean delays of less than approximately 
one billion years from star formation, then there could be on the order of only one SN la 
discovery at z > 1.4. However, a very meaningful test this survey could provide is if the 
TZiaiz) greatly increases with redshift, as it would if it were tied to the SFR//{( 2 ;), where one 
could expect at least few SNe la at z > 1.4 in this deep held survey. 

We have searched spans of images of the UDF, UDFP, and IRUDF observations for 
high redshift SNe and have found a total of four events, none of which were at z > 1.4. 
In §2 we describe the UDF+UDFP and IRUDF searches, show discovery and subsequent 
photometric data for the events found in the optical data, and present the evidence which 
show the relatively low redshift of these events. In §3 we compare the observed low observed 
yield at z > 1.4 to that which would be predicted from a few SN rate models. 


2. The Searches, Discovery Data, and ACS Photometry 

The UDF and UDFP images were searched by differencing stacks of images in the 
F850LP-band, assembled as they were collected over the duration of the deep survey. A 
total of nine multi-orbit epoch image stacks were created for the UDF held, and two were 
made for the UDFP. The hrst stack for each held was diherenced with an overlapping area 
in either the AGS GOODS mosaic in the case of the UDF, or with AGS images obtained 
for GO 9352 (A. Riess, PI) in the case of the UDFP. The subsequent stacks in both surveys 
were diherenced with their preceding stacks. Similarly, the IRUDF images were searched 
by diherencing two epochs of mosaic images, each 4 orbits in depth, in both the FllOIU 
and F160IU passbands. The two mosaic image stacks of the IRUDF were separated by 
approximately 77 days. In each the UDF, UDFP, and IRUDF searches, image stacks were 
produced using multidrizzle (Koekemoer et al. 2002) by using the median of several exposures. 
This successfully rejects cosmic rays without rejecting transients which last substantially 
longer than the length of one exposure (~ 1215 seconds). The UDF and UDFP were drizzled 
to 0."05 per pixel in these preliminary image stacks, and later re-drizzled to 0."03 per pixel 
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for the final release of the full-depth images. The IRUDF was drizzled to 0."09 per pixel, 
chosen to be a multiple of the hnal UDF pixel scale. 

Subtracted images were mined for candidate SNe by automated computer searches and 
close visual inspection of the images. The criteria for identifying potential SNe, and for 
rejecting possible confusion sources, were very similar to those used in S04. Each image 
stack of the UDF and UDFP would be sensitive to motions only larger than the point 
spread function (si O'M FWHM in F850LP) over the accumulation time of the image stack 
(roughly one week). The proper motion of any candidate would have to be less than 6 x 10“^ 
arcseconds hr“^ (1.4 x 10“^ degrees yr“^) to be misidentihed. This would exclude any 
foreground solar system object, as it would need to be farther than 247, 000 AU if orbiting 
the Sun at ~ 30 km s“^, which is well beyond the current estimates of the extent of the solar 
system. Slow-moving Galactic stars could also be generally excluded on similar hypothetical 
grounds. Variable stars are also easily excluded by their nature of being stars in all epochs 
of observation, and by generally appearing unassociated with galaxies in the held. Variable 
stars that, on occasion, do appear aligned with held galaxies are typically several orders of 
magnitude brighter than one would expect for SNe in these galaxies. 

What remained in the subtracted images are true extragalactic transients (e. g. super¬ 
novae and active galactic nuclei) and image artifacts such as misregistrations and variations 
in the point spread function. 

As the held is fairly dense with galaxies, and well resolved in the ACS images, there 
were no signihcant local misregistrations in the UDF and UDFP, except for stellar dihraction 
spikes which rolled with the change in orientation between epochs of the surveys. The point 
spread function (PSF) in each median combined image stack was largely stable and unahected 
by short-term, exposure-to-exposure changes in the PSF (e. g. telescope breathing or focus 
drift). Thus, we are conhdent that we have identihed only extragalactic transients from the 
UDF and UDFP subtraction images. 

Variable active galactic nuclei (AGN) were dehned as candidates that were within one 
pixel of their host nuclei. This is a rather conservative dehnition, as in principle, we would be 
capable of identifying transients with centroids a fraction of a pixel offset from the centroid 
of light in the nucleus of a galaxy. However, there were no detected candidates (down to the 
5cr threshold) that were < 0".05 from the centroid of the host galaxy the UDF and UDFP 
searches. Other investigations which probe for much smaller optical variations (1 to 3 (t) in 
galaxy nuclei do show evidence for low-level AGN in the UDF target held (Windhorst et ah, 
in preparation), none of which would be undeniably identihed with our diherencing method. 

The NIGMOS images, however, suhered from intrapixel sensitivity ehects. Gamera 3 
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is poorly sampled at ~ 0".2 per pixel, and can produce as much as a 30% flux variation 
depending if the peak of a PSF is centered on a given pixel, or is offset towards the edges 
of the pixel (Storrs et ah 1999; Lauer 1999; Hook & Fruchter 2000). This flux variation in 
the pixel response function (PRF) is not removed in the flat-held correction, and imposes 
an inherent limitation in these images. With many dithered exposures of the same held, 
it is possible to map the PRF (Storrs et ah 1999; Lauer 1999), and determine a correction 
to PSF photometry. Unfortunately, there are no easy means to deconvolve the PRF from 
the PSF. Therefore, differencing the two image stacks of the IRUDF necessarily resulted in 
a “checkerboard” pattern of under and over subtractions (measured at about ±15 — 20%). 
Faint and diffuse objects were generally cleanly subtracted. 

In Tables 1 and 2 we list the epoch stacks created for the UDF, UDFP, and IRUDF 
searches, with the mean date of the image stack, the total number of orbits, and the total 
exposure time of each stack. To asses the sensitivity and completeness of the UDF, UDFP, 
and IRUDF searches, we performed monte carlo tests with planted PSFs meant to represent 
false SNe. 

In the case of the IRUDF, a PSF was generated from a few stars in the images, and scaled 
to one count per second using zero points and aperture corrections provided by the NICMOS 
group at Space Telescope Science Institute. The measured aperture photometry was then 
corrected for intrapixel sensitivity variations following a prescription detailed in Storrs et al. 
(1999). 

We placed one false SN of a given magnitude on the centers of a randomly selected set 
of 50 detected objects (to > 5cr), subtracted the images, and then attempted to recover the 
false SNe by visual inspection (without prior knowledge of the locations of each planted false 
SN). This method was iterated to successively fainter magnitudes (in steps of 0.2 mag) until 
none of the planted SNe were recovered. The resulting histograms of percent recovered per 
magnitude bin is shown in Figure 1. 

The sensitivity in each epoch of the IRUDF was determined from the difference mag¬ 
nitude, Am, determined from the flux difference of objects in the residual frame, which is 
given by: 


mi = ZP- 2.5 X log(Fi) ^ Fi = 

m2 = ZP- 2.5 X log(F2) ^ ±2 = (1) 

Am = ZP — 2.5 X log(F 2 — F^i), 


where Fi and F 2 are the flux of the false SN in first and second image stacks, rrii and m 2 are 



the corresponding magnitudes, and ZP is the photometric zero point (total flux of 1 count 
per second). We use an analytical function to describe the efficiency in detecting objects of 
a given difference magnitude: 


T 

e{Am) — ^ g(Am-Amc)/S ’ 

where T is the maximum efficiency, Arric represents a cutoff magnitude where e{Am) drops 
below 50% of T, and S controls the shape of the roll-off. We And that T = 0.98, rric = 25.60, 
and S = 0.20 well describes the efficiency histogram for the FllOW passband, and T = 0.97, 
rric = 24.20, and S = 0.21 parameterizes the efficiency in the F160W passband (see Figure 1). 

In the F850LP passband of the UDF and UDFP surveys, false star tests show nearly 
identical efficiency function parameters as those used in S04, with T = 1 and S = 0.38, 
however with adjusted 50% efficiency cutoff magnitudes (me) corresponding to the 5cr sensi¬ 
tivity limits for the difference of a given pair of stacks (also shown in Tables 1 and 2). These 
limits are in good agreement with those expected from the exposure times of the search and 
template images using the ACS Exposure Time Calculators. The typical brightness thresh¬ 
old of detection for the F850LP template-search pairs was 27 mag. An illustration of the 
sensitivity threshold can be seen in Figure 2 where we have added six fake SNe (PSFs) at 
F850LP = 24, 25, 26, 27, 27.5, and 28 mag in a region of the 1224-1230 image stack. By dif¬ 
ferencing it with the 1212-1218 image stack, we can clearly detect the sources to a flux level 
of 27 mag, but fainter than 27 mag was difficult to identify without prior knowledge of where 
the fake SNe were. This is in good agreement with the estimated detection threshold (from 
the exposure times) of 27.1 mag for this template-search pair. The detection thresholds are 
bright enough to detect SNe la at peak at 2 ; < 2.2 in each passband. 

Despite the intrapixel sensitivity limitations of the IRDUF survey, interestingly, no can¬ 
didate SNe were discovered in the deep IRUDF imaging. However, four SNe were discovered 
in the UDF and UDFP images. For each, vega based aperture magnitudes were measured 
from difference images in the F606W, F775W, and F850LP bandpasses, using aperture 
corrections and photometric error estimations described in S04. In all cases, we used pho¬ 
tometric redshifts (phot- 2 ;) determined from the multi-wavelength GOODS data to estimate 
the redshifts of the host galaxies (Mobasher et al. 2004). Spectroscopic confirmation has not 
been obtained for these SNe. We also generally lack the photometric data and age constraints 
necessary to assuredly identify the SN types using the identification confidence scheme de¬ 
tailed in S04 and color selection methods described in Riess et al. (2004a). However, we can 
still use these techniques to reject combinations of SN type and redshift space, specifically 
we can reject the faint and red signature of SNe la at very high redshifts. The discovery 
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images for all four SNe are shown in Figure 3, and the positional and photometric data are 
listed in Tables 3 and 4. 

Three SNe were identihed in the area of the UDFP held. K0302-001^ was discovered in 
the GOODS follow-up images from GO 9352, and not detected seven months later (to within 
ha) in the UDFP images. The host of K0302-001 was very faint, with F6061U ~ 28 mag 
within a 0.2" radius, and very blue, virtually undetectable at F850LP > 27.5 in the same 
aperture. The galaxy could not be sufficiently detected in any of the deep multi-wavelength 
ground-based data sample assembled for the GOODS (spanning U though i7s-bands), and 
was only identihed in two AGS bands: F606W and F77h. The lack of photometric mea¬ 
surements of the host galaxy made it difficult to constrain its photometric redshift using the 
Bayesian Photometric Redshift (BPZ) method (Benitez 2000). The phot-;^ estimate derived 
of the host of K0302-001 lacks a signihcant peak and has a broad 95% conhdence interval of 
1.03 < < 2.22. However, the photometry of the SN is more illuminating. The magnitude 

and colors (specihcally the red F606hF — F850LP color of 3.33 mag; see Riess et ah 2004a) 
match those of the hve SNe la measured by Riess et ah (2004b) at or near maximum at 
2 : 1.3. Results from Riess et ah (2004a,b) demonstrate that ~ 95% of SNe with these 

photometric characteristics are correctly identihed as SNe la (with the other 5% being SNe 
Ic or Ib). A hrm conclusion we can make is that this SN is too bright and too blue to be a 
SN la (of any previously seen luminosity and color) at 2 : > 1.4. 

SN 20031t was discovered in the hrst UDFP stack (mean date 2003 Aug. 31), and 
undetected in the GOODS follow-up comparison images. It was also well detected in the 
second UDFP stack of mean date 2003 Sep. 12. The host galaxy was well detected in several 
passbands, and thus the phot- 2 ; for the host was well constrained at z = 1.0 (0.74 < ;2 < 1.26 
95% conhdence interval). We hnd that the photometry and colors of this SN were consistent 
with those for a SN la discovered ~ 80 days from maximum light at 2 : = 1.0, and inconsistent 
with tested SN la scenarios (varying age, light-curve shape, and extinction) at 2 ; > 1.4. 

SN 20031u was found in the second UDFP stack, and not detected in the hrst UDFP 
stack. The single F850LP-band measurement alone does not allow for a restriction in SN 
type and redshift space. However, the phot- 2 ; for the bright and well-detected host was 
;2 = 0.11 (0.0 - 0.25 95% conhdence interval), and therefore we can reject the possibility that 
this was a SN la at 2 ; > 1.4. 

The only SN detected in several epochs of the UDF target held observations was 
SN 2004R, discovered in the last 6-orbit stack obtained 2004 Jan. 13. It was not detected, to 


^This designation reflects the new International Astronomical Union (lAU) standard for possible faint 
supernovae. See http://cfa-www.harvard.edu/iau/CBAUSN.html 



within 5a, in the 12-orbit stack from 2003 Dec. 24 - 30. A review of the data obtained 2004 
Jan. 1-11 revealed that the SN was rapidly rising in the F850LP and F775hh-bands, but 
declining in the F606hh-band. The phot-z for the host galaxy was again well constrained at 
= 0.8 (0.56 < z < 1.04 95% conhdence interval). The light curves and colors were generally 
consistent with our Type Il-Plateau model, SN 1999em, caught prior to rest-frame 5—band 
maximum light, with a host extinction of ApsmLP ~ 1-5 mag (assuming a Galactic extinction 
law) at the photometric redshift. The colors and magnitudes of the SN are inconsistent with 
SN la models at z > 1.4. 

To summarize, our searching of the UDF, UDFP, and IRUDF resulted in four SNe of 
largely unknown types, but we can reject the possibility than any were SNe la at z > 1.4. 


3. Comparisons to Predicted Rate Models 

The sample of SNe la from the GOODS has shown for the hrst time a distinct rise in 
the SN la rate from 0.5 < z < 1.0 (Dahlen et ah 2004). These measurements are statistically 
inconsistent with a constant (or nearly constant) TZia{z), as inferred from measures of the 
SN rate at z < 0.5 (Blanc et al. 2004). 

An apparent discrepancy exists in the SN la rate measurements at z ^ 0.5. At this 
redshift, the SN la rate measured from the GOODS data are nearly twice as high as mea¬ 
surements made in the ground-based high- 2 ; supernova search programs (Pain et al. 2002; 
Tonry et al. 2003). The natural explanation would appear to be that completeness was not a 
primary goal for these ground-based programs, and was generally sacrihced for sake of favor¬ 
able SNe la with low background contamination for optimal spectroscopic conhrmation and 
precise distance measurements. Indeed, a recent, careful re-examination of the data from 
the Fall 1999 campaign of the High- 2 : Supernova Search project (Tonry et al. 2003), and the 
Fall 2001 IfA survey (Barris et al. 2004) shows preliminary evidence of many additional SNe 
resulting an increase by a factor of two in the SN la candidates, and thus a likely increase in 
the rates determined from these data to values consistent with the GOODS measurements 
(Brian J. Barris and John L. Tonry, in private communications). 

Perhaps the most intriguing result on the supernovae rates from the GOODS sample 
has been the dearth of SNe la discovered at 2 : > 1. Dahlen et al. (2004) have found that 
there is a steep decline in the rate at 1.2 < ;2 < 1.6 in comparison to the rate at z = 1.0, 
a result seemingly at odds with the lack of such a decline in the global star formation rate 
in this redshift range. However, S04 and Dahlen et al. (2004) have shown that TZia{z) is a 
reflection of the SFR[/( 2 :) which has been delayed by the assembly time function (or delay 
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time function), i. e. the time for SN la progenitors to go from formation to explosion. S04 
find that the difference in the evolution of the SFR[/( 2 ;) and the Tija{z) can be well modeled 
as a Gaussian delay time function with a mean delay of ~ 4 Gyr and a dispersion that is a 
small fraction of the delay. Subsequently Dahlen et ah (2004) hnd that based on the model 
delay, the fraction of stellar mass in the range of 3 to 8 solar masses which will explode as 
SNe la is about 5%. 

Past searches for SNe la in high-redshift clusters with HST by Gal-Yam et ah (2002) 
show rates that are similar to field SN la rates, indicating some similarity in their rate 
evolution. Unfortunately, these cluster rate measurement at 2 : 3> 0.1 are not precise enough, 
nor sufficiently well sampled in redshift space, to dehnitively compare their SN la rate history 
to those measured in normal (held) environments. 

An anecdotal anomoly in the consideration of the decline in TZia{z) at > 1 is the 
serendipitous discovery of SN 1997ff, found in only one re-observation of a pointing with 
HST and WFPG2 (Gilliland et al. 1999; Riess et al. 2001). The combined depth of the 
template-search pair for the SN 1997ff discovery image was only marginally deeper than for 
a single tile of template-search epoch of the GOODS. The SN itself, although likely magnihed 
by 0.3 mag (Benitez et al. 2002), was discovered well above the detection threshold for the 
tiny survey. The 5 month durration of the UDF-I-UDFP search and the factor of 2 increase in 
area of AGS over WFPG2 suggest that the number of such SNe la at 2 : > 1 found in the UDF 
should exceed those found by Gilliland et al. (1999) by a factor of ~5. Yet, when dealing 
with such small, pencil-beam surveys, the survey yields can be subject more to chance than 
expectations. The low number of expected SNe serve more as a probability for hnding SNe, 
governed by Poisson statistics. 

In the previous section we concluded that our searches of the UDF, UDFP, and IRUDF 
had failed to yield any SNe la at 2 ; > 1.4. To test the signihcance of this null result, we 
have simulated the yield expected in the multiple campaigns (image stacks) of the survey 
based on the best-£t 7lia{z) determined in S04. As in S04, the number of expected SNe per 
redshift interval was determined by: 


Nla{z) = nia{z%{z){l + Z)-^^AV{Z), (3) 

47r 

where tc is the “control time” probability function, which takes into account the survey 
efficiency with difference magnitude (described in §2), 0 is the area surveyed, and AV{z) 
is the volume element in a Az shell about 2 : for an assumed flat universe (Dfc = 0)^. The 


^The volume is V{z) = (47r/3)Z)p = (47r/3)D|/(l -I- z)^. This is the correct form of the equation, which 
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expected number of SNe la at 2 ; > 1.4 is then, 


2.6 

NUz > 1.4) = Nja{z)Az. (4) 

z=lA 

For comparison, we also calculated the yield expected from a 7lia{z) which is directly 
proportional (i. e. without delay) to the extinction corrected SFR; 7 ( 2 ;) model used in S04, 
but with an explosion efficiency for the SN la progenitor mass range of 10% (instead of 5%) 
as required to match SN la rate measurements at z < 0.2 by various authors (see Figure 4). 

As a straw man, we also computed the yield expected from a Tija{z) which proportional 
to the SFR/j:j( 2 ;) without any delay where we again assumed the progenitor efficiency to 
be 10% to match the low-z SN la rates. The expected number of SNe la from the tested 
models over all redshift ranges and at > 1.4 for each image stack of the UDF and UDFP 
(differenced with the prior image stack) are tabulated and summed in Tables 1 and 2. In 
Figures 5 and 6 we show the redshift distribution (total number in each redshift bin) expected 
by each Tlia{z) model, in the F850LP, FllOlF, and FlOOhF-bands. 

The best-fit TZiaiz) model predicts only 0.06 SNe la at > 1.4 in the UDF-I-UDFP, 
and 0.02 and 0.01 in the IRUDF FllOhF and F160hF-bands respectively, due to the long 
delay required by progenitors in the model. The universe was only 4.5 billion years old at 
2 : = 1.4 (assuming Ho = 70 km s“^ Mpc“^, him = 0.3, and IIa = 0.7), and few early forming 
progenitors would have had sufficient time to explode. We determined the significance of the 
yield given the model using the Poisson probability for observing a number of events. 




x\ 


(5) 


where x is the integer number of observed SNe la (x = 0 in this survey) and ^ is the expected 
number of events determined from the model. The absence of SN la events from this survey 
in the tested redshift range is an expected result for this best-fit TZja{z) model, with Poisson 
probabilities of 94.2%, 98.2%, and 99.1% assigned to this outcome in the F850LP, FllOhF, 
and F160hF-bands, respectively. However, a null result could also be expected from a TZia{z) 
proprotional to the SFR; 7 ( 2 ) (i.e., with no delay) with Poisson probabilities of 50.4%, 81.7%, 
and 86.7% for zero observed events in the F850LP, FllOlF, and F1601F-bands. Although 
the greatest expected yield would come from the model proportional to the SFR//{( 2 ;), a null 


has a typographical error in S04. 
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yield could still be expected in the tested redshift range with probabilities of 22.5%, 64.4%, 
74.8% in the respective passbands. 

It is interesting to consider what the relative agreement with the tested models would 
have been had a single SN la at 2 : > 1.4 been discovered in each search. If one high- 2 : SN la 
had been found in each the FllOhh and F160hh-bands, then it would have been very unlikely 
that S04 best-£t 7lia{z) represents the true SN la rate history, as such an outcome would 
be expected less than 2% of the time. The S04 best-fit Tlja{z) could therefore be rejected to 
> 98% confidence. However, one SN may have been expected in each IRUDF passbands for 
either the SFRt 7 ( 2 ;) or SFR/i:j( 2 :) models, although with Poisson expectation probabilities of 
12 — 28%, it would appear to be a low likelihood. Similarly, had one SN la at 2 : > 1.4 been 
discovered in the F850LP surveys, then S04 best-fit Tija{z) could be rejected at > 94%. 
But again the yield would be only marginally expected from either of the other two models, 
with a Poisson probability of ~ 35%. 

Although it appears as there was a statistical preference by the data for the best-fit 
7lia{z) model, it is clear that our null result cannot reject any of the tested models to a 
significant (greater than 99%) confidence. In fact, a model would need to predict five or 
more SNe la at 1.4 < ;2 < 2.4 to be rejected to greater than 99% confidence by the zero yield 
of this survey. Such a hypothetical model suggests that the average SN la rate in this high- 2 : 
range would be approximately 150 times larger than the SN la rate measured in the local 
{z < 0.1) universe. This is very inconsistent with trends inferred from measured SN la rates 
over any redshift interval, and is likely to be rejected based on theoretical considerations of 
the low [0/Fe] ratios it would predict for stars and the ISM of local galaxies. 

Such an increase in the SN la rate is plausible in the SN la rate models considered 
by Maoz & Gal-Yam (2004) in clusters of galaxies. Using e—folding delay time functions, 
described in detail in Madau et al. (1998) and Dahlen & Fransson (1999), with e—fold times of 
less than 1 Gyr, one could expect a SN la rate as high as 45 SNu in galaxy clusters at 2 ; ~ 1.5, 
following a burst of star formation in these clusters at ;2 = 2. However, to be consistent with 
the observed iron masses of galaxy clusters (assuming most of this mass is attributed to 
SNe la), the SN la rate history would likely need to drop to zero at 2 : < 0.7 — 1.0. This 
would be inconsistent with the non-negligible SN la rate measurements in galaxy clusters at 
0.04 < 2 < 0.08 (Reiss 2000). 

There is significant evidence for uncertainty in the estimates of galaxy volume densities 
for specific populations in the GOODS fields (Somerville et al. 2004), and the survey pre¬ 
sented in this paper covers a smaller area of the GOODS South field. However, the SFRt/( 2 ) 
model (in the 2 > 1 regime) used in this analysis was also determined from GOODS data, 
and therefore should be subject to the same variances. The SFR/ii;( 2 ) model, by contrast. 



was determined from several local surveys, and through theoretical modeling of the evolu¬ 
tion of the galaxy infrared luminosity function (Chary & Elbaz 2001). Therefore, we expect 
cosmic variance to distort the predicted SN la rate in this observed volume (from 15 — 20%) 
when applying the global average SFR/j:j( 2 ;) to this small held. 

It should be noted that our lower-bound on the redshift range tested, z > lA, was based 
on the redshift range where the completeness of the UDF-I-UDFP survey becomes signihcant 
over the GOODS SN survey. However, our results are not very sensitive to the precise 
dehnition of the uniquely high-redshift space of this survey. Additionally, other baselines 
and combinations of images in the UDF and UDFP surveys could have been selected which 
would have increased the expected number of SNe la in the desired redshift range. However, 
an increase by a factor of 5 or more would be necessary to begin to signihcantly rule out 
some of the tested models. Tests of varying cadences between visits and combining image 
stacks result in changes in the expected yield of the tested models of less than a factor of 
two in the considered redshift range. 

It appears that deepest-field-type, pencil beam surveys (i.e., the UDF or the Hubble 
Deep Field) with HST lack the volume necessary to accumulate a statistically signihcant 
sample to differentiate between the current plausible TZia{z) models in the redshift regimes 
in which they are uniquely sensitive. This is not surprising, considering that the increased ex¬ 
posure times over the GOODS did not signihcantly increase the luminosity distance probed. 
It would be necessary to repeat the UDF survey several times (or an increase of several 
times the survey area) in order to achieve expected yields for signihcant Poisson probability 
rejection of the tested models. This would not be economical with HST given the enormity 
of the UDF survey. 

Progress in determining the rates of SNe la and their associated timescales for assembly 
is most likely to come by “piggybacking” such studies on the now numerous and major ehorts 
to hnd and measure supernovae to constrain dark energy parameters. Such “cosmology- 
driven” SN surveys are designed as wide, open held surveys and are well-suited to deriving the 
supernova held rates. These surveys generally avoid galaxy clusters as search targets because 
half of the yield of SNe la found in such surveys are located behind the clusters (Gal-Yam 
et ah 2002; Reiss et al. 1998) where contamination by lensing is unavoidable. In addition, 
potential evolutionary diherences between cluster and held-born SNe la could add systematic 
errors to cosmological determinations. As a result, very deep rates and cluster rates of SNe 
la have been studied on a much smaller scale than held rates. 
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4. Summary 

The UDF, UDFP, and IRUDF observations yielded four SNe over the 5 month duration 
of the survey. We hnd the dearth of SNe la at high redshift in this survey to be consistent 
with the best-£t 7lja{z) from S04. However, it is also a likely result from a 7lja{z) which 
follows the SFR[/(; 2 ) at a higher progenitor explosion efficiency. More rapidly increasing 
T^iaiz), such as the Chary & Elbaz (2001) SFR/ij(; 2 ) model without delay, also cannot be 
significantly rejected by these data, because of the relatively small area probed in this survey. 
The general conclusion is that these data are well supported by the GOODS best-fit model, 
which shows that SNe la progenitors require substantial (approximately 4 Gyr) delays from 
star formation to produce events. However, other plausible SN la rate models cannot be 
excluded. Surveys such as the GOODS more than compensate for the loss in depth by the 
increased survey area, and are vastly superior programs for differentiating between SN la 
rate models. Future endeavors with HST such as the Cycle 12 programs by A. Riess and S. 
Perlmutter (Pis for GOs 9727 and 9728 respectively) and the Probing Acceleration Now with 
Supernovae (or PANS, GO-10189, A. Riess, PI) in Cycle 13 will allow for much more precise 
measure of the SN la rate in the z > 1 regime. In combination with the GOODS data, they 
will provide the best available measures of the empirical distributions of SN la delay times 
and the comparison to star formation rate models. 

We thank Tomas Dahlen for his beneficial discussions. We also thank Steven Beckwith, 
Roger Thompson, Bahram Mobasher, Anton Koekemoer, Louis Bergeron, Rogier Windhorst, 
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Table 1. Predicted yield of SNe la from UDF+UDFP 


Field 

Mean Date 

# of 

Exp. Time 

5cr lim. 

Baseline 


TZia{z) models 




Orbits 

(seconds) 

(mag) 

(days) 

S04 best-fit 

SFRc/ 

SFRir 

GOODS 

2003 Feb 11 

5 

10500 






0924-1002 

2003 Sep 28 

18 

43740 

27.3 

229 

1.013 (0.017) 

1.256 (0.206) 

3.425 (0.444) 

1002-1008 

2003 Oct 05 

18 

43740 

27.4 

7 

0.031 (0.002) 

0.065 (0.023) 

0.169 (0.050) 

1010-1014 

2003 Oct 12 

18 

43740 

27.4 

7 

0.033 (0.003) 

0.076 (0.032) 

0.195 (0.071) 

1016-1029 

2003 Oct 22 

20 

48600 

27.5 

11 

0.052 (0.005) 

0.129 (0.059) 

0.326 (0.128^ote. — The ex] 

1204-1211 

2003 Dec 07 

16 

38480 

27.3 

46 

0.190 (0.014) 

0.429 (0.173) 

1.097 (0.®B(J with the expec 

1212-1218 

2003 Dec 15 

18 

41940 

27.4 

8 

0.035 (0.002) 

0.075 (0.027) 

0.195 (0.Cfia)re survey. 

1224-1230 

2003 Dec 27 

12 

29160 

27.1 

12 

0.049 (0.004) 

0.105 (0.038) 

0.272 (0.084) 

0101-0111 

2004 Jan 06 

8 

19440 

26.9 

10 

0.040 (0.002) 

0.076 (0.022) 

0.202 (0.049) 

0113-0114 

2004 Jan 13 

6 

14580 

26.7 

8 

0.032 (0.002) 

0.057 (0.015) 

0.153 (0.033) 






UDF = 

1.475 (0.051) 

2.268 (0.595) 

6.034 h-290) 

GOODS 

2003 Feb 11 

2 

4200 






0831-0901 

2003 Aug 31 

2 

4600 

26.1 

202 

0.595 (0.006) 

0.651 (0.056) 

1.766 (0.125) 

0911-0913 

2003 Sep 12 

9 

20700 

26.8 

12 

0.058 (0.003) 

0.114 (0.034) 

0.302 (0.075) 






UDFP = 

0.653 (0.009) 

0.765 (0.090) 

2.068 (0.200) 






TOTAL = 

2.128 (0.060) 

3.033 (0.685) 

8.100 (1.490) 


Table 2. Predicted yield of SNe la from NICMOS UDF 


Field 

Mean Date 

# of 
Orbits 

Exp. Time 
(seconds) 

5cr lim. 
(mag) 

Baseline 

(days) 

S04 best-fit 

F-iaiz) models 
SFRu 

SFRir 

FllOW 

Stack 1 

2003 Sep 07 

4 

10752 






Stack 2 

2003 Nov 23 

4 

10752 

25.6 

77 

0.341 (0.017) 

0.605 (0.170) 

1.598 (0.374) 

Fimw 
Stack 1 

2003 Sep 07 

4 

10752 






Stack 2 

2003 Nov 23 

4 

10752 

24.2 

77 

0.234 (0.012) 

0.481 (0.216) 

1.166 (0.429) 





o 


(3) Aou0!O!jj3 


Fig. 1.— Efficiency of the IRUDF survey in recovering false SN placed in the centers of 
galaxies in the held. The fraction of recovered fake SNe per magnitude bin is shown as 
a histogram in the FllOlF and F160fF passbands, and is approximated by the function 
e{Am) oc (1 + 
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F850LP=24.0 S/N = 100 


0 

F850LP=27.0 S/N = 7 


Fig. 2.— Example of the average survey sensitivity in the F850LP-band. A region of the 
1224-1230 image stack is shown after differencing it with the previous image stack. Fake 
SNe added at various magnitudes. Below the image are radial prohles (solid points, lines 
are Moffat hts to the prohles) of the 24 mag and 27 mag fake SNe, with approximate S/N 
of 100 and 7 respectively. SNe above ~ 27.1 mag (which corresponds to the 5a limit) were 
virtually undetectable in this survey. 






19 


2003 Fob 11 


2003 Aug 31 ^ 


K0303>001 


N 


. i 


I 


1 


2003 Aug 31 


2003 F«b 11 


8M 200311 


2003 Sup 12 


2003 Aug 31 


8N 30031U 


10 * 


2.5" 


2004 Jan 13 

2003 Sap 28 

8N 2004R 


► 


- 

5" 


Fig. 3.— Discovery images for SNe K0302-001, 20031t, 20031u, and 2004R (first column), 
shown with template images without the supernovae (middle column), and subtracted resid¬ 
ual frame (last column). North is up and East to the left in all images. 






Fig. 4.— The Type la supernova rate history, as measured by several authors (points). 
Vertical bars represent statistical errors, and horizontal bars on Dahlen et al. (2004) points 
represent redshift bin size. Also shown are the tested models for the TZia{z): The best-£t 
delayed SFR model (solid line), a SFR-like model (dashed line), and Chary & Elbaz (2001) 
SFR model (dotted line). 
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Table 3. Positional Data 


SN 

R. A. (2000) 

Deck (2000) 

N (arcsec) 

E (arcsec) 

Redshift 

K0302-001 

03:32:37.10 

-27:56:53.6 

-0.27 

0.13 

1.3 (1.2 - 1.4) 

20031t 

03:32:42.88 

-27:55:52.5 

0.45 

0.20 

1.0 (0.74 - 1.26) 

20031u 

03:32:36.17 

-27:55:01.4 

-0.28 

0.43 

0.11 (0.0-0.24) 

2004R 

03:32:41.30 

-27:46:13.6 

-0.10 

-0.11 

0.80 (0.56 - 1.04) 


Note. — Offsets are from the nucleus of the host galaxy to the SN. Maximum likeli¬ 
hood photometric redshifts are listed with the 95% conhdence intervals in parentheses. 
Coordinates supersede those announced in the lAU Circular. 







Table 4. Photometric Data 


SN 

Filter 

JD+ 2,450,000 

Magnitude 

K0302-001 

F850LP 

2680.22 

24.38 (0.03) 


F775W 

2680.96 

25.20 (0.03) 


F606W 

2680.25 

27.71 (0.12) 

20031t 

F850LP 

2890.58 

25.85 (0.06) 



2895.71 

26.19 (0.04) 


F775W 

2885.32 

27.26 (0.11) 


F606W 

2887.52 

30.2 (0.8) 

20031u 

F850LP 

2895.71 

26.33 (0.04) 

2004R 

F850LP 

3006.95 

27.54 (0.26) 



3014.59 

26.79 (0.13) 



3015.49 

26.87 (0.10) 



3017.13 

26.49 (0.06) 


F775W 

3005.89 

27.91 (0.26) 



3009.21 

27.09 (0.12) 



3013.70 

27.48 (0.07) 



3017.45 

27.12 (0.05) 


F606hF 

3000.40 

27.82 (0.07) 



3007.75 

27.85 (0.06) 



3016.28 

28.04 (0.08) 



3019.15 

28.20 (0.09) 


Note. — Magnitudes are given in the Vega-based 
system and are listed with their photometric errors in 
parentheses. 








